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The effect of journal eccentricity on the static and dynamic performance 
of a water lubricated, 5-recess hybrid bearing is presented in detail. The 
hydrostatic bearing has been designed to operate at a high speed and with a 
large level of external pressurization. The operating conditions determine the 
flow in the bearing to be highly turbulent and strongly dominated by fluid 
inertia effects. The analysis covers the spectrum of journal center 
displacements directed towards the middle of a recess and towards the mid-land 
portion between 2 consecutive recesses. Predicted dynamic force coefficients are 
uniform for small to moderate eccentricities. For large journal center 
displacements, fluid cavitation and recess position determine large changes in 
the bearing dynamic performance. The effect of fluid inertia force coefficients 
on the threshold speed of instability and whirl ratio of a single mass flexible 
rotor is discussed. 


INTRODUCTION 


Hydrostatic and hybrid(combined hydrostatic and hydrodynamic) journal 
bearings (HJBs) have recently been the subject of increased attention because of 
their potential applications as support elements in cryogenic high speed 
turbomachinery. These bearings have been proposed because of their high radial 
stiffness, accuracy of positioning and good vibration damping characteristics. 
Despite these attractive features, the hydrodynamic action at high speeds may 
bring the potential of half speed whirl instability. Thus, it is of importance 
to determine the limits of safe operation for such bearing applications 

To provide the levels of load support needed in high speed machinery 
operating with cryogenic fluids requires large external pressure supplies. 

These unique operating conditions give rise to large film Reynolds numbers, and 
cause the flow within the bearing film to be highly turbulent and strongly 
dominated by fluid inertia effects both at the recess edges and in the thin film 
land regions. Conventional analytical treatments are confined to the numerical 
solution of the Reynolds equation for turbulent flow, and thus, can not 
interpret correctly the actual complex phenomena occur ing in HJBs. 

*The present research effort has been sponsored by Rocketdyne, by the NASA 
Center for Space Technology, and by the State of Texas Advanced Technology 
Program. 
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NOMENCLATURE 


b Recess circumferential length [m] 

c radial clearance [m] 

C XX’ C XY’ C YX’ C YY damping coefficients [Ns/m] 


D 

d 


F , F 
x’ y 


fj = 0.001375 


f B = 0.001375 


Orifice discharge coefficient 
2 R . Bearing diameter [m]. 

Orifice diameter [m] 

Fluid film forces in X and Y directions [N]. 
f 5xl0 5 \ 1/2>65 i 


1 + 


1 + 


R, 


5x10 


Rr 


5 .1/2.65, 


Turbulent friction factor at journal 
surface 


Turbulent friction factor at bearing 
surface 


H/c = l+e^cos© + e y sin0 Dimensionless film thickness. 


I 1 ,I 2 ,I 4 


Recess depth [m] . 

Inertia coefficients for stability analysis. 

K XX ’ k XY ’ K YX ’ ^YY Stiffness coefficients [N/m] 

Equivalent bearing stiffness coefficient [N/m]. 
Half of rotor stiffness [N/m]. 


K 


eq 


K 


k , k 
x* y 


k J’ k B 


L 

1 

N 


rec 


(l/2)(kj+k.g) dimensionless shear parameters in 
circumferential and axial directions 

fjRj, fgRg* Turbulent shear parameters at journal and 
bearing surfaces. 

Bearing axial width [ m ] . 

Recess axial length [m]. 

Number of recesses on bearing. 


M XX’^XY’ M YX ,M YY Inert ^ a or added mass coefficients [Kg]. 


M 


Half of rotor mass [Kg]. 


Fluid pressure [N/m 4 ]. 
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2 

External supply pressure [N/m ] . 

f H (V.n) dr.. Recess flow into film lands [m3/s] 

pS2cR/y Nominal circumferential flow Reynolds number. 

3 2 

pc (P -P )/y R Pressure flow Reynolds number. 

S 3 . 

(Re^) (c/R) Modified pressure flow Reynolds number. 

2 2 1/2 2 2 1/2 
( p/y)H[ (U-S2R) +V Z ] / , ( p/y)H[U+V*] A/ . 

Reynolds numbers at journal and bearing surfaces. 

Fluid film resistance torque [Nm] 

mean flow velocities in circumferential and axial dirs. 

3 

(H j .+H)A i .+V s Total recess volume [m ] . 

3 

Volume of orifice supply line [m ]. 

Inertial coordinate system. 

relative effective fluid inertia mass at support. 

2 

Fluid compressibility factor [m /N]. 

(e ,e )/c. Dimensionless journal eccentricities 
7 in X & Y directions 

2 2 

pS2R /c (Ps-Pa). Rotational speed parameter. 

00 / 52 . Whirl frequency ratio, 
s s 

Whirl frequency ratio without fluid inertia effects. 
Excitation or whirling frequency f 1/sec] 

Whirl frequency at instability threshold [1/sec]. 
Rotational speed of journal [1/s] 

Threshold speed of instability [1/s]. 


Redecliff and Vohr (1969), Heller (1974), and Deguerce et al. (1975) have 
presented finite difference solutions to the turbulent flow Reynolds equation. 
Reported numerical flow rates and hydrodynamic forces are lower than 
experimental measurements but considered satisfactory. Ar tiles, Walowit, and 
Shapiro (1982) presented a well documented numerical method for prediction of 
the static and dynamic performance characteristics of turbulent HJBs. The 
analysis includes the local effect of fluid inertia at the recess edges, but the 
fluid compressibility effect at the recess volumes is neglected. The turbulence 
model employed by Artiles seems to give inconsistant results in the transition 
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from pressure dominated flow to pure shear flow at high speed operation. 

The influence of fluid inertia at the recess edges and hydrodynamic 
pressure generation within the recess for high rotational speeds have been shown 
and discussed in detail in the experimental work of Chaomleffel and Nicholas 
(1986). The magnitude of the measured edge pressure drop appears to be dependent 
on the journal rotational speed and the recess depth. Comparisons of the 
experiments with predictions based on a finite difference solution to the 
turbulent Reynolds equation with fluid inertia effects are also included in the 
same reference. The numerical calculations for load and flow rate were found to 
be always greater than the measurements, and predicted edge pressure drops were 
determined to be smaller than those observed. 

SanAndres( 1989a) has developed a full inertial numerical analysis for 
accurate prediction of the performance characteristics of orifice compensated 
turbulent HJBs. The analysis departs from the mean flow Navier-Stokes equations 
and includes the effect of recess volume fluid compressibility. The numerical 
results presented by SanAndres (1989b) indicate that fluid inertia at the film 
lands acts as an additional flow resistance, and thus, causes reduced flow rates 
and greater hydrodynamic effects when compared to an inertialess solution. 
Inertia force coefficients of large magnitude in bearings handling large density 
fluids and with low values of the clearance ratio (c/R) are reported. The full 
inertial analysis developed by the author shows a considerable improvement over 
conventional analyses based on the classical Reynolds equation. 

The present paper introduces further numerical predictions for a high speed 
HJB currently under testing. The effect of journal eccentricity on the static 
and dynamic performance of the sample bearing are throughly detailed. A simple 
analysis to determine the influence of fluid inertia on the threshold speed of 
instability for a simple flexible rotor supported on hybrid bearings is 
presented. 


ANALYSIS 


Figure 1 presents the geometry of a hydrostatic journal bearing. Fluid 
passes across a restrictor orifice and a recess area of relatively large 
clearance, and then flows through the film lands to the bearing discharge. The 
performance characteristics of the bearing are governed by momentum and 
continuity considerations in the fluid film lands as well as flow continuity 
through the bearing recesses. The former defines the pressure and velocity 
fields over the bearing lands, whereas the latter provides values for the 
pressure in the recesses. 

San Andres (1989a) presented a derivation of the appropriate equations for 
fully turbulent flow in hybrid bearings. The analysis considers the effect of 
fluid inertia (advection and temporal terms) at the film lands and a local 
Bernoulli type inertia effect at the recess edges. For an isoviscous 
incompressible fluid, the continuity and momentum equations for the turbulent 
bulk flow on the film lands are given as 
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The turbulent shear parameters k and k are taken as local functions of 
the friction factors relative to X the journal and bearing surfaces (Hirs, 1973, 
Nelson, 1987). 


The continuity equation at the recess is defined by the global balance 
between the flow through the orifice restrictor, the recess outflow into the 
film lands, and the temporal change of fluid mass in the recess volume. The 
recess continuity equation is then given by 


nd 


2(P -P ) 
v s r' 


1/2 


3V 3P 

r r 

= Q + — — + V 6 — 

at at 


(4) 


where (3 represents the fluid compressibility factor at the recess volume. The 
recess outflow into the film lands, Q , is determined by 

integration of the mean flow velocity vector along a closed path T encompassing 
the recess boundary. 

For hybrid operation a pressure rise occurs within the downstream portion 
of the recess (Chaomleffel et al.,1985). This region is considered as a 
one-dimensional step bearing. The local acceleration of fluid from relatively 
stagnant conditions at the recess to a high velocity at the film lands causes a 
sudden pressure drop at the recess boundary. The pressure at the entrance to 
the film lands is modeled by simple Bernoulli type relations based on the theory 
developed by Constantinescu (1975). 

San Andres (1989a) presented a numerical procedure for effective solution 
of the nonlinear governing equations (1) through (5). The numerical method is 
based on the SIMPLEC algorithm of Van Doormal & Raithby (1984) for iterative 
solution of the film land equations, and it is combined with a Newton-Raphson 
algorithm for estimation of recess pressures and satisfaction of the recess flow 
constraints. The flow field is divided into a zero-th and first order fields 
describing the equilibrium condition and perturbed dynamic motion, respectively. 
Solution of the zero-th order field gives the steady state bearing 
characteristics such as flow, load and friction torque. The first order field 
determines the dynamic force coefficients for the excitation frequency of 
interest. 
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RESULTS 


A Hydrostatic-Hybrid Bearing Test program is under development at Texas 
A&M University. The purpose of the program is to develop an experimental data 
base and analytical capabilities for turbulent flow hydrostatic bearings. The 
objectives of the program are to design and build a test apparatus, 
systematically test the bearings, and develop an analysis to predict the bearing 
characteristics and anchor it to the test data. 

The hydrostatic bearing test apparatus has been designed and is currently 
under preliminary testing. Table 1 presents a description of the 5 recess test 
hydrostatic bearing considered as a test case. The working fluid is hot water at 
a temperature of 71.1C(160F), and the maximum rated inlet supply fluid pressure 
is 7.23 MPa(1050 psi). The rotor speed ranges from 0 to 500Hz(30000 rpm), and 
the bearing and journal surfaces are considered perfectly smooth. For the 
conditions shown in Table 1, the flow within the bearing is highly turbulent 
with dominance of hydrodynamic effects at the highest rotational speed. The 
large value of the modified Reynolds number, Re , shows that fluid inertia at 
film lands will have a considerable influence on^the test bearing static and 
dynamic performance characteristics. 

The orifices of the test bearing have been sized for a concentric pressure 
ratio equal to 0.5 at the maximum operating speed of 30000 rpm. This condition 
implies that the predicted pressure drop between the supply and recess pressures 
is equal to the pressure drop between the recess and the outside of the bearing. 
The selected concentric pressure ratio is not necessarily an optimum for the 
present bearing configuration but was chosen on the grounds of common practice 
for establishing bearing stiffnesses which are not too far from optimum. For an 
orifice discharge coefficient equal to one, the present full inertial analysis 
predicted an orifice diameter equal to 1.496mm. 

A complete analysis has been carried out to determine the effects of 
journal center eccentricity on the static and dynamic performance 
characteristics of the test bearing at 30000 rpm and an external supply pressure 
equal to 6.54 Mpa(950psi). The journal center displacements ranged from the 
concentric position to eccentricity radii equal to 90% of the film clearance in 
the directions toward the middle of a recess and toward the mid-point in the 
lands between two consecutive recesses. Figure 2 shows the predicted bearing 
flow rate and the film friction torque for the range of journal displacements 
considered. The flow rate is a maximum at the concentric position and decreases 
as the eccentricity increases. The friction torque is a minimum at the centered 
position and increases rapidly as the journal center displacement approaches the 
bearing clearance. The increment in torque and reduction in flow rate are due to 
increased flow resistance for growing eccentricities. The results show the flow 
and the torque to appear closely as even functions of the eccentricity. The 
maximum reduction in flow rate at the largest eccentricity considered, 

(e/c)=0.9, is approximately equal to 22% of the flow at the centered position, 
while the increment in torque is about 30%. 

Figure 3 shows the calculated recess pressures for different values of the 
static journal eccentricity. Negative eccentricities, e , denote journal 
displacements to the midland between recesses 1 and 5 located at 36° and 324°, 
respectively. Positive eccentricities^ e , denote journal displacements toward 
the middle of recess 3 located at 180°. Figure 4 shows the entire pressure field 
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over the journal surface for selected values of the journal center displacement. 
At the concentric position the recess pressures are identical and the bearing 
has null load capacity. The pressure field for null eccentricity is depicted in 
Figure 4(c). 

For displacements toward the middle of a recess, Figure 3 shows that 
maximum recess pressures are attained at recesses 2 and 3 located 
in the converging portion of the film thickness. The large pressures attained at 
the second recess are caused by the effect of journal rotation and will generate 
a substantial film force normal to the journal displacement. For eccentricity 
radii larger than 0.5, cavitation was found to appear in the land region around 
recess 4 where the film thicknesses is diverging. Figures 4(d) and 4(f) show the 
pressure field for eccentricities e x equal to 0.4 and 0.7, respectively; note 
that for large journal displacements there is a significant hydrodynamic effect 
between recesses 2 and 3. 

For displacements toward the midland between recesses 1 and 5, Figure 3 
shows that maximum recess pressures are located at recess 5 in the convergent 
portion of the film thickness while low recess pressures are attained at 
recesses 2 and 3 in the divergent film thickness region. Fluid cavitation is 
predicted to occur for eccentricities larger than 0.6. Figures 4(a) and 4(b) 
show the pressure field for eccentricities equal to -0.7 and -0.4, 
respectively. Note from Figure 3 that for large journal displacements toward the 
land between the recesses, the pressure at recess 5 exceeds the supply pressure 
indicating that there is backflow through the orifice. These results are to be 
taken with caution, since the analysis presumes the flow-orifice relation to be 
valid for any flow direction. However, experimental evidence presented by 
Redecliff and Vohr (1969) shows that, for high eccentricity operation, recess 
pressures larger than the external supply pressure are indeed possible. 

Figure 5 shows the predicted film forces, -F and F , parallel and normal 
to the static journal displacement, respectively. X For th^ journal speed 
considered, 500 Hz(30000 rpm) , the hydrodynamic effects are dominant and 
determine a normal force F larger than the hydrostatic force F . The angle 
between the resultant load^and the eccentricity vector is equal X to 56° and 
remains relatively constant for the entire range of eccentricity radii 
considered. For small to moderate eccentricities, the film forces are linear 
with the static displacement and will determine constant stiffness coefficients. 
For large static displacements toward the middle of a recess, the curvature of 
the film forces is upward due to fluid cavitation. However, for displacements 
toward the mid-land, the film forces curvature is first inward reducing the film 
forces and then, outwards increasing the forces. This seemingly strange effect 
is explained in terms of the pressure field not being able to develop 
substantially in the convergent portion of the film thickness due to the 
presence of recess 5, and thus pushing water back-flow through the orifice. 

Figures 6 through 8 show the calculated synchronous stiffness, damping, and 
inertia coefficients for the hydrostatic bearing operating at 30000 rpm. The 
rotordynamic coefficients are uniform for small to moderate journal eccentricity 
radii. The cross-coupled stiffness coefficients shown in Figure 6 are 
substantially larger than the direct coefficients and demonstrate the dominance 
of hydrodynamic effects on the bearing dynamic performance. For large static 
journal displacements, the bearing direct stiffness (K ) increases dramatically 
due mainly to the effect of fluid film cavitation. Figures 7 and 8 show the 
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direct damping and inertia coefficients to be one order of magnitude larger than 
the cross coupled damping and inertia coefficients. The magnitude of the direct 
inertia coefficients is large due to the large density of the fluid and the low 
clearance to bearing radius ratio. At the concentric position, the direct 
inertia coefficient is equal to 17.9 Kg and determines for synchronous 
excitation at 50^ Hz(30000 rpm) a reduction of almost 88.3% in the dynamic 
stiffness (K -toM ). This result proves that fluid inertia is of utmost 
importance in the dynamic performance of hydrostatic bearings handling fluids of 
large density. 

The calculated rotordynamic coefficients for an excitation frequency equal 
to the rotating speed are not affected by the compressibility of the fluid at 
the recess volume due to the low recess depth to clearance ratio of the designed 
test bearing. However, for supersynchronous excitations well above the rated 
journal rotational speed, compressibility effects become important and may cause 
significant losses in damping and increased values of direct stiffness. 

The results obtained for the test bearing dynamic coefficients prompted the 
question of the test rotor-bearing stability characteristics. Appendix A 
presents the criteria derived for the threshold speed of instability and whirl 
frequency ratio of a simple flexible rotor supported on identical bearings and 
includes the effect of the inertia force coefficients. This effect is 
significant in the present application. 

Figure 9 presents the calculated whirl frequency ratio ( <|>) versus the 
journal eccentricity for the rotor-bearing system studied. The calculated 
results show that the whirl ratio is independent of the rotor flexibility and 
equal to 0.492 for small to moderate static journal center displacements. For 
the largest positive eccentricities (journal displacement towards the middle of 
a recess), the whirl ratio drops sharply to a value equal to 0.36. However, for 
large displacements toward the middle of a land (negative eccentricity e ) the 
whirl ratio increases well above 0.5. This abnormal effect is due to the loss in 
direct stiffness and the inability of the bearing to attain high pressures at 
the film lands due to orifice backflow. The results presented clearly show that 
the dynamic characteristics of the test hybrid bearing are no better than for a 
plain journal bearing. 

Figure 10 presents the threshold speed of instability for the test 
rotor-bearing system. The rotor mass is equal to 30.84 Kg and the shaft is 
considered as rigid. The figure includes the predicted onset speed £ if fluid 
inertia effects are neglected. A net reduction in critical speed of almost 25% 
is observed to occur due to the large influence of the inertia force 
coefficients on the stability response of the test system. 


CONCLUSIONS 

The effect of journal eccentricity on the static and dynamic performance 
of a water lubricated, 5-recess hybrid bearing is presented in detail. The 
hydrostatic bearing has been designed to operate at a speed of 30Krpm and a 
large level of external pressurization(6.64MPa) . These operating conditions 
determine the flow in the bearing to be highly turbulent and strongly dominated 
by fluid inertia effects. 
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The numerical predictions are based on the full turbulence-inertial 
analysis of SanAndres( 1989a) . The analysis covers the spectrum of journal 
center displacements directed towards the middle of a recess and towards the 
mid-land portion between 2 consecutive recesses. Bearing flow rate decreases 
and friction torque increases as the eccentricity increases. Film forces are 
linear with the eccentricity for small to moderate static journal center 
displacements. For large eccentricities, film forces are affected by fluid 
cavitation and present some unusual features. For large static displacements 
towards the mid-land it is observed that back flow through the orifices occurs. 
This unexpected behavior degrades considerably the static as well as the dynamic 
performance of the bearing. 

The dynamic force coefficients for the test bearing are predicted to be 
constant for small to moderate values of the journal eccentricity. However, for 
large static displacements toward the mid-land it is observed that the direct 
force coefficients are substantially reduced. Direct inertia force coefficients 
of substantial magnitude are predicted in the present application. The levels 
obtained are a direct consequence of the large density of the fluid considered. 

Hydrostatic bearings operating in hybrid mode at high rotational speeds 
present similar stability characteristics as plain journal bearings. For a 
single mass flexible rotor supported on HJBs, the threshold speed of instability 
is reduced by the rotor flexibility and the fluid inertia effect. The 
cross-coupled inertia coefficients are shown to degrade the stability response 
of the bearing by reducing the whirl frequency ratio. 


APPENDIX A 

Linear Stability Analysis 

Criteria for the threshold speed of instability and whirl frequency 
ratio of simple rotor-bearing systems are well known in the literature (Lund, 
1965, Chen et al., 1985). These analyses usually consider a rigid or flexible, 
single mass rotor supported in bearings which possess stiffness and damping 
characteristics. The influence of inertia force coefficients in the stability of 
these rotor-bearing systems has not been accounted due presumably to the 
general assumption that these coefficients are negligible in most circumstances. 
However, the calculated inertia force coefficients for the present application 
are rather large compared to the mass of the test rotor, and thus, their effect 
on the stability of the system is regarded as important. 

Consider a single rotor of mass 2M supported on a pair of identical 
bearings by a massless shaft of stiffness 5K r . Instability sets in at a 
rotational speed when the journal starts a forward precessional motion at the 
whirl frequency w = <j> Q . From the solution of the homogeneous system of 
equations for the motion of the rotor and bearing supports, it is determined 
that the whirl speed w and the whirl frequency ratio 4> are easily obtained from 
the following quartic equations: 
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If fluid inertia effects are eliminated, then the coefficients a., I., 

I 2 and 1^ are equal to zero, and the stability formulae derived are identical to 
the criteria presented by Lund (1965) and Chen(1985). From equation (A.l) it is 
determined that the fundamental effect of fluid inertia is to reduce the 
threshold speed of instability if the ratio of masses 0 ^ is large. 

For operation at the concentric position (e =c =0), the direct dynamic 
coefficients are identical and the cross-coupled coefficients differ in sign. At 
this condition, the whirl frequency ratio is obtained from the expression? 

4> 2 S2 2 M +<f)SJC - K =0 (A. 5) 

r s xy T s xx xy 


where it is noted that the whirl ratio is reduced if the cross-coupled inertia 
coefficient Mxy is negative. 
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Table 1. Description of 5 recess hydrostatic vater bearing under 

development for TAMU Hybrid Bearing Research program (1989) 
Bearing geometry and operating conditions. 
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Figure 2. Predicted flow rate and friction torque for test hydrostatic 
water bearing as a function of static journal eccentricity. 
Pressure supply equal to 6.54MPa(950psi) and journal speed 
equal to 500Hz( 30000 rpn). Orifice dianeter equal to 1.496 mo. 



Figure 3. Predicted dimensionless recess pressure for test hydrostatic 
water ©earing as a function of static journal eccentricity. 
Pressure supply equal tc 6.54HPa(950psi) and journal speed 
equal to 5Q0Hz( 30000 rca). Orifice diaaeter equal to 1.496 ma 
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Eccentricity (e /c ) 

5. Film forces parallel(-F ) and normalfF ) to journal static 

displacements for test hydrostatic vat^r bearing. Ps-6.54 MPa, 
30000 rpa, orifice diaaeter-1.496 mm. 



Figure 6. Predicted stiffness coefficients for test hydrostatic water 

bearing versus journal static eccentricity. Ps-6.54 MPa, 30000 
rpa, orifice diaaeter-1.496 aa. 
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Eccentricity Ex 


Figure 9. Whirl frequency ratio + for test rotor-hybrid bearing system 
versus journal static eccentricity. 



Figure 10. Threshold speed of instability 2 for test rotor-hybrid 
bearing systea versus journal static eccentricity. 
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